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ABSTRACT 



Recently, several new correlations between gamma ray burst (GRB) observ- 
ables have been discovered. Like previously well established correlations, they 
challenge GRB models. Here we show that in the cannonball (CB) model of 
GRBs, the newly discovered correlations have the same simple kinematic origin 
as those discovered earlier. They all result from the strong dependence of the 
observed radiations on the Lorentz and Doppler factors of the jet of highly rel- 
ativistic plasmoids (CBs) that produces the observed radiations by interaction 
with the medium through which it propagates. 

Subject headings: gamma rays: bursts 
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Introduction 



Despite the enormous complexity and diversity of gamma ray bursts (GRBs) and 
their afterglows, various well established correlations between GRB observables were found 
during the past years. Such correlations challenge the theoretical models of GRBs. Most 
of these correlations have been neither predicted nor explained by the standard fireball 
model of GRBs, which has been extensively employed to explain the GRB phenome non (for 



reviews see, e.g., 



Meszaros 2002 



Zhang fc Meszaros 2004 



Piran 2004 : 



Zhang 20071 ). In the 



cannonball (CB) model of GRBs, these correlations were shown to be a simple consequence 
of the strong dependence of GRB observables on the Lorentz factor 7 and Doppler factor 5 
of the highly relativistic jet of plasmoids (CBs) whose interaction with the medium along its 
path produces the observed radiations (Dar & De Riijula 2000,2004; Dado et al. 2007,2009). 

Recently, many new correlations were discovered between pairs of observables 
characterizing th e prompt gamma r ay and early optical emissions in a sample of GRBs rich 



in X-ray flashes (ILiang et al. 20101 ) . Here we show that in the CB model all these new 
correlations also follow from the strong dependence of the GRB observables on the Lorentz 
factor 7 and the Doppler factor 5 of the highly relativistic jet of CBs. 

In particular, from their selected GRB sample, Liang et al. (2010) inferred tight pair 
correlations between the initial Lorentz factor To of the relativistic ejecta in GRBs and 
other GRB observables, such as the isotropic equivalent gamma ray energy Ei Son of the 
prompt emission, the peak luminosity L Py o of the early optical emission, and the peak time 
t' of the optical emission in the GRB rest frame. The fireball model itself does not explain 
the physical origin of these newly discovered correlations, nor the origin of other well 
established correlations between GRB observables. Moreover, in the fireball model, the tight 
r — E is0tl correlation yields a peak energy E va that is p ractically independent of Ej iSiLa . 



in contradiction with the well established Amati relation (lAmati et al. 2002 



Amati 20061 ) . 
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However, the initial value of the Lorentz factor of the r elativistic ejecta t hat produced the 



GRBs was not truly measured. Instead, it was inferred ( ILiang et al. 20101 ) from an assumed 
relation, T oc lEj ^/jt') 3 ] 1 ^ 8 . This relation follows from the standard fireball model 
expression for T (ISari fc Piran 19991 ) . neglecting its weak dependence on the circumburst 
density and on the conversion efficiency of the relativistic kinetic energy of the ejecta 
to radiation. The validity of this fireball model relation has never been proven. Thus, 
hereafter, we shall distinguish between the true value of the bulk motion Lorentz factor 7 
of the relativistic ejecta and r , its alleged value inferred from fireball modeling of GRB 
data. We shall show that, in the CB model, the strong dependence of E isoa and t' and 
consequently of T on the true initial values of the Lorentz factor 7 and the Doppler factor 
5 also yields the tight pair correlations between the fireball model parameter Tq and the 
GRB observables that were discovered by Liang et al. (2010) and could not be explained 
by the fireball model. 



2. Kinematic origin of correlations in the CB model 



In the cannonball (CB) model of GRBs (Dado et al. 2002,2004,2009), gamma ray 
bursts and their afterglows are produced by the interaction of bipolar jets of highly 
relativistic plasmoids (CBs) of ordinary m atter with the radiation and matter along their 



trajectory (jShaviv &: Par 1995 



Par 19981 ). Such jetted CBs are presumably ejected in 



accretion episode s on the newly f ormed compact ste l lar object in core-collap se supernova 



(SN) explosions (IPar et al. 19921 ; |Par fc 



compact objects in close binary syst ems ( 



phase transitions in compact stars (IPar 1998 



laea 19991; IPar fc Pe Ruiula 2000T). in me rger of 



Goodman et al. 1987: 



Shaviv fc Par 1995) and in 



Par & Pe Ruiula 2000 



Pado et al. 20091 ). 



GRBs are not standard candles because of the diversity of their central engines and 
environments. But, because of the large bulk motion Lorentz factor 7 of the jet of CBs, 
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their emitted radiation at redshift z, which is observed at a small angle 9 relative to 
the direction of the jet, is boosted by a large Doppler factor 5 = I/7 (1 — (3 cosO) and 
collimated through relativistic beaming by a factor 5 2 . Moreover, the time difference dt in 
the observer frame between the arrival of photons emitted by the point-like CBs at two 
different points along their path, which are separated by a distance dr' and time dt' = dr'/c 
in the progenitor's rest frame (hereafter a prime indicates an observable in the progenitor's 
rest frame), is shortened (aberrated) according to dt — dt'(l + z)/ r y5. The large Doppler 
boosting, relativistic collimation and time aberration produce correlations between GRB 
observables, despite their dependence on the CBs' intrinsic (rest frame) properties and on 
the environment along their trajectories (which produce a significant spread around these 
simple kinematic correlations). 

The redshift z of the GRB location is measurable, and the dependence of the GRB 
observables on redshift can be taken into account explicitly, unlike their dependence on the 
values of the Lorentz factor and the viewing angle of the jet, which can only be inferred 
with model-dependent assumptions. However, the strong dependence on 7 and 5 can be 
used to correlate triplets of independent observables without knowing the values of 7 and 
5. Moreover, several observables depend on the same combination of 7 and S that result in 
pair correlations. Finally, due to selection effects in the observations, various observables 
depend strongly only on 7 or 5, which also yields pair-correlations. In particular, the 
dependence on viewing angle of the jet can be eliminated in two general cases: For 7 2 1 
and small viewing angles 9 2 <C 1, the Doppler factor satisfies <5~2 7/ (I+7 2 9 2 ) to an excellent 
approximation. For # 2 7 2 <1 the Doppler factor decreases rather slowly with increasing 
9. But, for 6> 2 7 2 3> 1 the Doppler factor decreases with increasing viewing angle like 9~ 2 , 
and the observed fluence of gamma rays, which in the CB model is amplified by a factor 
5 2 due to relativistic beaming, decreases like [1 + 7 2 6> 2 ]~ 2 . The geometrical probability to 
view a bipolar GRB from a small angle 9 increases like (1 — cos9) m9 2 /2 and the product 
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9 [1 + j 9 ] has a maximum when ^9 =1. Con sequently, 5 = 7 for the most probable 



viewing angle 9 = l/j of GRBs (jShaviv h Par 19951 ). 



In 'soft' GRBs with j 2 9 2 > 1, such as X-ray flashes (XRFs) with 7 2 # 2 >10, the 
dependence on the exact value of the Lorentz factor can be ignored compared to the strong 
dependence on the Doppler factor. This yields a slightly different correlation. Thus, we 
shall derive the CB model correlations for the above two different situations, for 5 ~ 7 and 
for r y 2 9 2 S> 1, i.e. 5 <C7 where the dependence on 7 can be neglected compared to that 
on 5. For a mixed population of ordinary GRBs, soft GRBs and XRFs one may expect 
an approximate correlation obtained by using the average index of the two power-law 
correlations. (Alternatively, one can derive triple correlations, i.e., correlations that involve 
three independent observables that each of which depends both on 7 and 9). Below we 
derive the the CB model pair correlations that correspond to those discovered empirically 
by Liang et al. (2010) and we compare the two in the text and in Table 1. 



3. Pair Correlations 

3.1. Correlations between temporal parameters 

In the CB model, the prompt optical and afterglow emissions of GRBs are emissions 
from the jet while it crosses the wind/ejecta blown from the progenitor into the interstellar 
medium (ISM). Until the CB sweeps in a relativistic mass/energy comparable to its rest 
mass, 7 and 5 stay put at their initial values. Because of time-aberration, all time measures 
of the prompt emission pulses/flares, such as its rise-time t r from half maximum to peak 
value, the peak-time t p after the beginning of the flare/pulse, the decay-time td from peak 
value to half maximum, and its full width at half maximum (FWHM) t w , are proportional 
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to (1+z)/j5. Hence, 

logU = a>ij + logtj , where i,j = r, d,p, w. ( 

This explains the origin of the 'universal' power-law index ~ 1 of the 'power-law 
correlations' (proportionality) between the temporal parameters of the prompt 7-ray 
pulses. Moreover, the CB model also predicted the pulse shape and the ratios between 
the temporal parameters of the prompt G RB pulses. They yield the simple relations 



(bar fc De Ruiula 2004 



Dado et al. 2009h : 



td~2.39t r , i.e., logtd~0.38 + logt r , 
t(i~1.41tp, i.e., logtd~0.15 + logt p , 
t r ^0.59t p , i.e., logt r ^—0.23+logt p , 
t w ^2.00t p , i.e., Iogt w tt0.30 + logt p , 
t w ~3.39t r , i.e., log t w tt0.53+logt r , 
t w &lA2t d , i.e., logt w ^0.15 + logt d , 

where time is measured in seconds. These relations are well satis fied within observatio nal 
errors by well resolved GRB peaks/flares as was found, e.g., in (jKocevski et al. 20031 ) 
through empirical parametrization of these peaks/flares. In the CB model these correlations 
are valid approximately also for the prompt optical flares as indeed was found by Liang et 
al (2010): 

logt d = (0.48 ± 0.13) + (1.06 ± 0.06) logt r , 
logt d = (-0.09 ± 0.29) + (1.17 ± 0.11) logt p , 
logt r = (-0.54 ± 0.22) + (1.11 ± 0.08) logt p , 
logt w = (0.05 ± 0.27) + (1.16 ± 0.10) logt p , 
logt w = (0.61 ± 0.11) + (1.05 ± 0.05) logt r , 
logt w = (0.15 ± 0.02) + (0.98 ± 0.01) logt d . 

Note, however, that Liang et al. (2010) have selected for their analysis GRBs with 
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a single optical peak that could be described by the empirical formula of Kocevski et al. 
(2003) for GRB pulses. However, in most GRBs, the prompt optical emission is probably a 
sum of unresolved ove rlapping; fla res, as suggested by s everal bri ght GRBs such as 50820A 



071031 ( 



(iVestrand et al. 2006^. 080319B flWozniak et al. 20091 ). 061007 flRvkoff et al. 2009h and 



Kruhler et al. 20091 ) . where the large photon statistics allowed a good temporal 



resolution with rapid-response ground-based telescopes. In nearby soft GRBs (GRBs with 
a relatively small E' p ) and XRFs, the individual flares are wider and observed later because 
of a much smaller Doppler factor, while the duration of the central engine activity remains 
the same, independent of viewing angle. Consequently, in far off-axis GRBs, such as XRFs, 
the individual flares strongly overlap and produce an effective single peak. Such peaks are 
hardly res olved into overlapp ing flares even with large rapid-response telescopes such as 



GROND flDado fc Par 2009al ). Thus, the sample of GRBs in Liang et al. (2010) consists 
mainly of soft GRBs and XRFs with a single optical peak at a relatively large peak-time, 
which allowed their measurement with rapid response telescopes that are larger but slower 
than the robotic telescopes. 

Note also that because the peak energy E Pt l of the prompt 7-ray pu l ses and of the earl v- 
time X-ray flares is proportional to r y5/(l + z) (IDar fc De Rujula 2004 iDado et al. 20091 ). 
all the above temporal parameters of the prompt 7-ray, X-ray and optical peaks/flares are 
also inversely proportional to E p of the 7-ray emission: t r ocl/E p ^, t p (xl/E Pjl , td(xl/E p ^, 
t d (xl/E pa . 



3.2. Pair correlations during the optical emission 



In the CB model, the optical e mission is domina ted by synchrotron radiation (SR) and 



its lightcurve is given by (see, e.g. 



Dado et al. 2009l .and references therein), 



F„[t] ocn^/^V 3 " 1 ^, 



(2) 
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where n[t] is the density along the CB trajectory, R is the CB radius and /3 is the spectral 
index of the SR. Typically, the spectral index in the optical band is 0o ~ 0.5 at early 
time and gradually approaches 0x ~ 1 at late time. The peak luminosity of the prompt 
optical emiss ion flare is obtained when the CB reaches the peak density of the progenitor's 
wind/ejecta ( iDado fc Par 2009al ) while 7 and 5 stay put at their initial values. Thus, for 
the most probable viewing angle of GRBs #^1/7, i.e., £^7, the CB model predicts for the 
peak optical luminosity L Pj o OC7 4 , and t' w oc I/7 5 ~ I/7 2 for the FWHM in the GRB rest 
fra me. These depend encies yield the correlation L Pt o(x[t' w ]~ 2 , in agreement with that found 



in ( ILiang et al. 20101 ) and reported in their equation (11) and in Table 1 below. Similar 



correlations are expected between L pQ and t r , t d and t p , respectively . 

Note, however, that in the CB model, XRFs that are far off-axis GRBs have initially 
7 2 # 2 3> 1. The deceleration of the CBs in XRFs yields S(t), which first rises slowly 
as j(t) decreases with increasing time until it reaches a maximum when j(t)9 = 1, 
i.e., when 5(t) = r y(t). Hence, for a constant density ISM, a typical early-time optical 
spectral index f3 Q — 0.5, and S(t) ~ y(t), Eq. 02]) yields L p Q oc [7(t)] 4 . Hence, for 
t' w oc l/y(t) 5(t) ~ [7(^)]~ 2 , the CB model predicts the correlation L Pj0 oc [C]~ 2 - 



3.3. Correlations involving total energy, peak energy and peak luminosity 



In the CB model the peak energy E p ~ and the isotropic equ ivalent gamma ray ener gy 



Eiso,-y of the prompt emission peaks and early-time flares satisfy (IDar &: De Rujula 20001 ) 



[1 + z) E pa oc 7 5, E ison oc 5 3 . 



(3) 



For 5^7, Eq. (EJ) yields (1+z) Ep^ocy 2 , £'j SOi7 oc7 3 and the correlations (1+z) E pa oc Ef^. 
In soft GRBs and XR Fs, i.e., GRBs with a la rge viewing angle, the CB model predicts 



1/3 

l + z) E Pt ry oc 5 oc E^ (IDar fc De Rujula 20001 ) . For a mixed population of soft and hard 
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GRBs, it yields the mean correlation 

(1 + z) E vn oc E\Hfl\ (4) 



consistent with the Amati correlation (lAmati et al. 20021 ; lAmati 20061 ). which was 
discovered empirically. 

In the CB model, where the early optical emission is dominated by synchrotron 
radiation with a canonical spectral index flo ~ 0.5, Eq. ([2]) predicts that ordinary GRBs 
with 5 p» 7 have a peak optical luminosity L p Q oc 7 4 and consequently L p Q ocEf^. In 
soft GRBs and XRFs where the dependence on 7 at early time can be neglected compared 
to the dependence on 5, the resulting correlation is L p>0 oc 5 7 ! 2 oc E 7 Jl r Thus, the mean 
correlation in a mixed population of soft and hard GRBs that is predicted by the CB model 



is. 



L P ,o oc £J 0)7 , (5) 

which is in agreement with the power-law correlation with an index 1.17 ± 0.13 found by 
Liang et al. (2010) and reported in their equation (12) and in Table 1 below. 

The equivalent isotropic optical energy of the prompt optical flare that dominates 
E is0:0 is roughly given by L p t' w (x^~ 1 ^ 2 S 5 ^ 2 . Thus, in ordinary GRBs where S « 7, the CB 
model predicts E is0) o OC7 2 , i.e., E is0) o ocEf^. In soft GRBs, E is0t ooc5 5 ^ 2 oc [E iso ^] 5 ^ 6 , and 
the mean effective correlation that follows is, 

E lso ,o oc <f /12 , (6) 

which is in good agreement with the correlation found by Liang et al. (2010) and reported 
in their equation (14), and in Table 1 below . 



3.4. Triple correlations 



Many correlations between a triplet of GRB observables can be derived using their 
strong dependence on 7 and 5. For instance, the peak equivalent isotropic luminosity of 
prompt emission 7-ray pulses and early-time X-ray flares satisfies, 

-^p,7 ~ Ei so ^ /t w OC Eiso^y Ep ,y . (7) 

In the CB m odel, the break time of a canonical X-ray afterglow satisfies 
( IDado et al. 20091 ). tb,x ^ (1 + z )/($ 2 7)- Then, using the relations in Eq. (j3J), the 
following correlation is obtained, 

^ a l/(£ p [4o l7 ] 1/3 ). (8) 

The three observables in the triple correlations depend usually also on other jet parameters 
and on the burst environment. These dependences can produce a significant spread around 
the predicted correlations. 



4. Origin of the tight T -E isoa correlation 

Liang et al. (2010) reported a tight correlation, T ~ 182 [E iso ^/10 52 erg] 25 , between 
r of the ejecta and E ison . However, T was not the measured initial Lorent factor but 
rather a value inferred from the relation 

r ~i92[^ so , 52i7 /(g 3 ] 1/8 ; (9) 

where £^0,52,7 = Ej so>1 /(1 0^_erg} and t' r is in seconds. This relation follows from a fireball 



model expression for T (ISari fc Piran 19991 ) after neglecting its weak dependence on the 



unknown circumburst density and on the fraction of the relativistic kinetic energy of the 
ejecta that is converted to radiation. 
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Substitution of the CB model relations E isoa oc 5 3 and if oc 1/(7 5) into Eq. (E} yields 
T = 7 3 / 8 5 6 / 8 . Consequently, for a mixed population of ordinary GRBs and XRFs, the CB 
model predicts the following correlations compared to the observed (obs) ones: 
T cc [^ SO , 7 ]°- 31±0 - 075 obs: T oc [£ JSO , 7 ]°' 25±0 - 02 
T oc [L Pi0 ]°' 25±0 ' 03 obs: T oc [L PiO ] a20±a03 

T OC [ty-0-66±0.09 obs . ro K ^]-0.59±0.05_ 

Note that for a GRB population dominated by XRFs and soft GRBs, where the dependence 
on 7 can be neglected compared to that on 5, the CB model predicts: T oc [E iso ^] - 25 , 
T oc [L Pi0 ] - 21 , and T oc [t' p }-°^. 



5. Summary and conclusions 

In a long series of publications we have demonstrated that the cannonball model 
of GRBs predicted correctly the main observed properties of GRBs, including the well 
established correlations between GRB observables, and can reproduce successfully their 
broad band lightcurves of GRBs and their afterglows from o nset until very late tim es, 



despite their eno r mous complexity and diversity (see, e. 



Dado et al. 2009 



Dado k Par 2009a 



Par fc De Riijula 2004 



Pado fc Par 2009bl and references therein) . In this 



paper we have shown that in the CB model all the newly discovered pair correlations 
between gamma ray burst observables, reported in Liang et al. (2010) as the previously well 
established correlations, are a simple consequence of the strong dependence of the GRB 
observables on the Lorentz factor and viewing angle of the highly relativistic and narrowly 
collimated jets whose interaction with the medium along their path produces the observed 
radiations. 

One of the new discovered correlations which was reported in Liang et al. (2010) was 
the surprisingly tight relation, To — 182 [i?j SOj7 /10 52 erg] 0,25 . The authors pointed out in their 
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paper that "there is no straight-forward theory that predicts this relationship between To 
and Ej sn J\ and that this tight correlation is inconsistent with the well established Amati 
relation ( lAmati 20061 ) if the prompt gamma ray emission in GRBs is produced by the 
internal shock mechanism of the fireball model. 

The tight correlation between T and E is0tl , however, was based on values of T inferred 
in Liang et al. (2010) from a fireball model relation r oc [E is0jJ / '(£p) 3 ] 1//8 and measured 
values of Ei soa and t' p , and not from reliable measurements of the initial Lorentz factor 7. 
Using the CB model dependences of E isoa and t' p on 7 and 5 one reproduces the power-law 
correlation between T and E ison with the observed index 0.25 for a population rich in soft 
GRBs and XRFs. Indeed, the GRB sample that was used in Liang et al. (2010) to infer 
the tight correlation contains a large fraction of XRFs (e.g., 060904B, 070318, 070419A, 
071010A, 080330, 080710, 070208) and soft GRBs. We conclude that the tight correlation 
between Tq and Ei Son is that expected in the CB mod el, as well as the Amati r elation (Eq. HJ), 



which actually was predicted by the CB model (see 



before it was discovered empirically ( lAmati et al. 2002 



Par fc De Ruiula 2000 



Eq. (40)) long 



Amati 20061 ). In the case of the 



standard fireball model, where the prompt emission pulses are produced by synchrotron 
emission from internal shocks, the observed tight correlation between Tq and Ei SOil yields 
essentiall y E pn that is con stant for different E isoa values, in contradiction with the Amati 
relation (ILiang et al. 20101 ). This, perhaps, is not surprising in view of the fact that the 
standard fireball model was never shown to predict correctly the shape of the prompt 
emission gamma ray pulses and their spectral evolution, nor to explain convincingly their 
typical photon energy and total energy. 

Finally, we would like to caution that the new correlations discovered by Liang et al. 
were inferred from a sample of selected GRBs: Only GRBs with a single optical peak that 
could be modeled with the empirical formula of Kocevski et al. (2003) seem to be included 
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in the sample. Bright GRBs with clear early time multipeak optical emission were not 
included in the GRB sample. A single late-time peak in fact could be a late-time flare or 
the blended sum of unresolved peaks. Moreover, soft GRBs and XRFs usually have a slow 
rebrightening of their optical and X-ray afterglows that probably has a completely different 
origin - an initial rise in S(t) as j(t) decreases due to deceleration. Thus, most XRFs and 



soft GRBs show after a prompt emission flare (s) a slowly rising afterg 



power-law decay similar to that of ordinary GRBs (IDado fc Par 2009d j 



ow with a late-time 
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Table 1. Comparison between the power-l aw index k of the 



power-law correlations 



log[Oi\=ciij + klog[Oj] which was reported in (ILiang et al. 20101 ) for various pairs (Oi,Oj) 

of GRB observables and its value predicted by the CB model. t r , t p , and t w stand, 
respectively, for the rise-time, peak-time, decay-time, and total width of the optical peak. 



Correlated Pair 


k [obs] 


k [CB model] 


{td-i tr) 


1.06±0.06 


1.00 


(td, tp) 


1.17±0.11 


1.00 


(tr j tp) 


1.11±0.08 


1.00 


{twi tp) 


1.16±0.10 


1.00 


(t-wi tr) 


1.05±0.05 


1.00 


(tw j td) 


0.98±0.01 


1.00 


(Lp,o,t p ) 


-2.49±0.39 


-2.00 


{L Pj o,f w ) 


-2.00±0.32 


-2.00 


(-Lp^Oy Ei S0) ry) 


1.17±0.13 


1.25±0.08 


(tp, Ei so ^) 


-0.40±0.07 


-0.50±0.17 




0.74±0.10 


0.75±0.08 




0.51±0.06 


0.50±0.17 


(roj Ei Son ) 


0.25 ±0.03 


0.31 ±0.07 


(fo) Lpp) 


0.20 ±0.03 


0.25 ±0.04 


(r ,g 


-0.59 ±0.05 


-0.66 ±0.09 
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